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I. INTRODUCTION 



Brownian motion is one of the paradigms of statistical mechanics. Following the seminal 
works of Einstein, Langevin and Smoluchowsky, a detailed mathematical analysis can be 
made on the basis of Langevin or Fokker-Planck equations. While Brownian motion and 
its "time derivative" , Gaussian white noise, are stochastic processes of fundamental impor- 
tance, the dichotomous Markov process (see, e.g., has its own virtues and interest. First 
- and this is our central contention here - systems driven by dichotomous noise can often 
be described in full analytic detail. Second, dichotomous noise reduces to white shot noise 
and Gaussian white noise in the appropriate limits Third, it can either mimic the effects 
of finite correlation time of the noise, or it may directly provide a good representation of 
an actual physical situation (such as, for example, thermal transitions between two config- 
urations or states). Finally, it has the advantage that it can easily be implemented as an 
external noise with finite support. 

Most of the results for dynamics driven by dichotomous Markov noise are limited to 
systems with a single scalar variable x{t). The rate of change of this variable switches at 
random between the "+" dynamics, x(t) = f+{x), and the "-" dynamics, x(t) = f-{x) (the 
dot stands for the temporal derivative). This dynamics can be described by the following 
stochastic differential equation: 

i^t) = i±^^^I^ + i±^^^l^I^at), (1) 

where ^(t) is a realization of the dichotomous Markov process taking the values ±1, with 
transition rates between these values equal to and k_ respectively. Several physical ap- 
plications, including persistent random walks, nuclear magnetic resonance, chromatography, 
Taylor dispersion, and kinetic theory, have been discussed when the speeds f+ and /_ are 
constants, and exact time- dependent results can be obtained in this case [3|. Somewhat 
surprisingly, the full time-dependent solution for linear /+(x) and f-{x) is not available, 
see for a detailed discussion. In the case of nonlinear dynamics, steady state properties 
can be calculated and have notably been studied in detail in the context of noise induced 
.....ion. B and .o,se ..duced phase ..an.tion. g Specific dy.a,n,c properties Kave 
also been obtained. A significant effort has gone into the calculation of first passage time 
moments [? I and of transition rates, cf . Kramers' rate for thermal escape and resonant 
activation j9|. Furthermore, when f+{x) and f-{x) are periodic another important dynamic 



property, namely, the asymptotic drift velocity, can be extracted from steady state results 
by working with periodic boundary conditions. This technique is of particular interest in 
the context of Josephson junctions and Brownian motors lOl lUl, llSf • However, contrary 
to the claims made in some of these papers, the problem of the first passage time moments, 
and the related issue of finding the asymptotic drift velocity was not solved in the most 
general case. Indeed, with a few exceptions ll^ 13, all the results that have been 
obtained exclude the cases when one or both of the "±" dynamics have unstable fixed points, 
and thus do not consider the possibility of crossing unstable fixed points in the long time 
dynamics. The technical subtleties were first highlighted and discussed in detail in a recent 

n 

paper [19|| . In view of the broad applicability and importance of dichotomous Markov noise, 
a comprehensive review of the results for the asymptotic drift velocity in periodic systems 
driven by such a noise is called for and is the subject of this paper. As a first illustration 
we also apply these results to the calculation of the drift for the rocking ratchet. 

II. SOLVING THE MASTER EQUATION 

The calculation of the asymptotic drift velocity in periodic systems is most easily carried 
out by starting from the master equation, equivalent to (Q), for the probability densities 
P+(x,t) and P_(a;, t) to be at x at time t, if ^ = +1 and —1 respectively: 

dP (r f) f) 

— ^ = -— [Uix)p+ix,t)] - KP4x,t) + k^P4x, ,t), 

dP (x t) d 

' ' = [f-ix)P-ix,t)] - k^P4x,t) + k+P4x,t). (2) 

We assume that f+{x) and f-{x) are continuous functions of their argument and that they 
are both periodic, i.e., 

f+{x) = f+{x + L) and = f^{x + L), Vx. (3) 

In order to extract the long-time average drift speed, it is sufficient to study the steady state 
properties of Eq. for x G [0, L] with periodic boundary conditions. To show this, we 
introduce the steady state quantities P{x) = P+(x)+P_(a;) (which represents the probability 



density for being at x regardless the value of ^) [27[, and p{x) = k+P^{x) — k^P„{x). From 
the summation of the two equations in Q one immediately concludes that for the asymptotic 
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(steady) state the probability flux J associated with P(x), 

J = ^_±LM±^JMp^,) + (4) 

is a constant. The asymptotic drift velocity is then simply given by: 

{x) = [ [Uix)P+{x) + f-ix)P^{x)] dx = LJ. (5) 

By substracting the equations in Q (multiplied, respectively, by fc+ and one obtains, 
in the asymptotic steady state, the following first-order differential equation for p{x): 

^ {Kk.[U{x) - f-{x)]P{x) + [fc+/+(x) + k.f.{x)]p{x)} + (A;+ + k.fp{x) = 0. (6) 

Equations (@)) and (jHI) have to be solved by imposing the conditions of continuity for P{x) 
and p{x) (or, at least, the condition that they have no more than integrable singularities) on 
[0,L], periodicity, P{x) = P{x + L) and p{x) = p{x + L),Wx, and normalization for P{x), 

/ P{x)dx = 1. 
Jo 

All these elements allow the determination of the unique steady state solution P{x) , 
the stationary probability flux J, and hence the corresponding dynamic quantity of interest 
to us, namely, the asymptotic drift velocity (x). As was pointed out in detail in [l^, the 
situation is entirely different, both physically and mathematically, when one or both of 
the "±" dynamics have unstable fixed points, depending on whether the system can cross 
(or not) these unstable fixed points in the long time limit. In fact, while this issue was also 
mentioned in j3] and , a full and detailed discussion was first given for a specific example 
with multiplicative dichotomous noise in 2. In the following we focus on the presentation 
and discussion of the final results for the asymptotic probability density and drift velocity, 
relegating the technical details to the Appendix. 



III. MAIN RESULTS 



Since we assume that the functions f±{x) are continuous, and in view of the periodicity, 
fixed points in the separate dynamics x = f+{x) and x = f-{x) will always appear in pairs. 
In the following we will present the final results for the three simplest cases that can occur, 
namely no fixed points, cf. Section IIII Al one of the dynamics has two fixed points and the 
other none. Section IIIIBl and both dynamics have two fixed points. Section IIII CI 
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A. No fixed points 



The stationary probability density P{x) reads: 

(x) [U{x)-f.{x)] 



P{x) 



LUix)f4x) 

fX+L 



X 



X exp 



dz 



exp 
k+ + k 







dz{k+IU{z)+k.lf.{z))]-l 



Mz)-nz) \2[Uiz) - f4z)] 



dw{k+/U{w) + k./Uw)) 



X 



and the mean asymptotic velocity is given by: 



(x) 



L 



X 



exp 

x+L 



dz{k+l!Az) + k-lf-{z)) 



f+[x)f-{x) 



dz 



k+ + k. 



Uiz) + Uz) 
U{z)-Uz) ' \2[U{z) - f4z)] 



+ 



(7) 



X exp ( - / dwik+/Uiw) + k./f4w))j I . (8) 

Here (...)' stands for the derivative with respect to x. Note that the above expressions for 
both P{x) and (x) display the required symmetry between the "+" and "-" dynamics. The 
above expression reduces to the one given earlier in the literature for the particular case of 

n 

additive dichotomous noise [f+{x) — /_(x)=constant, cf. 



B. One of the dynamics has two fixed points in [0,L) 



We suppose that the "+" dynamics has two fixed points xi < X2 in [0, L), xi being stable 
= 0, < 0], while X2 is unstable [f+{x2) = 0, /+(a;2) > 0]. The "-" dynamics 

has no fixed points. 

As explained in detail in Appendix IVIBI and in fl9^ , there exists exactly one solution 
P{x) of the Eqs. ^ and (jH)) that is physically and mathematically acceptable, and for 
X e [a;i,a;i + L] it is given by the following expression: 



Pix) ^ f 



k+ + k^ 



dz sgn 



X2 



[U{z)-f.{z) 



U{z)-Uz) \2[U{z)-Uz)] 
X exp ( - / dw{k+/U{w) + k^/f^{w)) ] . 



(9) 



As discussed in Appendix IVIBl P{x) is finite and continuous throughout the interval 
{xi, Xl + L), with, in particular, at the unstable fixed point X2 

hm P{x) = hni P{x) = — — — — — — — — ^— -r. (10) 

^/^2 x\x2 L /_(x2) [A;+//|(x2) + 1] 

Note that positivity of P{x) implies that f-{x) and (i) must have the same sign. Hence the 
direction of the asymptotic drift is the same as that of the "-" dynamics, which is physically 
obvious. 

At the stable fixed point Xi, P{x) is either continuous when k^/\f'j^{xi)\ > 1, i.e., 
r P( ^ r p( ^ (fc+ + fc_)/|/;(xi)| - 1 

hm P{x) = hm P{x) = — — r-j- — r-rf- — -r, (11) 

x\x, x/{x,+L) L /.(xi) [/c+/|/;(xi)| - 1] ^ 

or divergent but integrable when A;_|_/|/_^(xi)| ^ 1 [namely, k^/\f'j^{xi)\ < 1 corresponds to 
a power-law divergence, while A;+/|/^(a;i)| = 1 corresponds to a "marginal" logarithmic-like 
integrable divergence]. This result is consistent with the physical intuition that probability 
density builds up near a stable fixed point, especially when the switching rate is low. 

The normalization condition / P(x)dx = 1, together with Eq. (0) lead to the 

Jxi 

following expression for the mean velocity: 



xi+L 



dx 



Xl 



f+{x)f-{x) 



dz sgn 



X2 



U{z)f.{z) 



X 



k+ + k. 



U{z) + f-{z) 
U{z)-f4z) ' \2[U{z)-f4z)] 



+ 



X 



X exp 



dw{k+/U{w) + k./f.{w)) 



This expression reduces to the one mentioned in 
additive dichotomous noise. 



(12) 



for the particular case of a symmetric 



C. Each of the alternating dynamics has two fixed points in [0,L) 

Consider as above xi < X2 as the stable, respectively unstable fixed points of the "+" 
dynamics, and suppose now that the "-" dynamics also has two fixed points, x^ < x^. 
Depending on the relative positions of these four points, as well as on their nature (stable 
or unstable), two different types of situations might occur. 

(i) In all the situations in which xi has as direct neighbor another stable fixed point, it 
is physically clear that the asymptotic dynamics settles into a random alternating motion 
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between these points, so that they dehmit the interval in which the steady-state probabihty 
density is non-zero [5|. Obviously in this case (i) = 0, i.e., there is no net flux of the 
particles, while P{x) is given by 

' ^ dz{K/U{z) + kJUiz))) , (13) 



P x) 



for X and xq lying in between the two stable flxed points, while the constant C is determined 
through the normalization condition for P{x). 

(ii) The cases that correspond to an alternation of the stable and unstable flxed points 
are of more interest to us, as they lead to a nontrivial behavior of P{x) and to a nonzero 
flow of the particles. Without loss of generality we suppose that < Xi (s) < x^ (u) < X4 
(s) < X2 (u) < L, with the following results: 



P(x) 



L 



U{x)-f4x) 



dz sgn 



U{z) + f.{z) 



X 



X 



Uiz)-f_iz) V2[/+W-/-(^)]. 
X exp dw{k+/f+{w) + k_/f.{w)) 

for X E (xi, X4) 



(14) 



(x) 
L 



f+{x) - f-{x) 



k, + A;_ 



dz sgn 



X2 



U{z) + f.{z) 



X 



U{z)-f4z) V2[/+(^)-/-(^)]. 

X exp dw{k+/f+{w) + k-/f.{w)) 

for X G (x4, Xi + L) 

P(x) is continuous throughout (xi,X4) and (X4, xi + L) and, in particular, at the unstable 
flxed points X2 and X3 it takes, respectively, the values (jlOj) and 

The periodicity of P(x) is connected with the behavior at the stable flxed point xi that 
was already discussed in the previous section, see Eq. (fTTj). 

Concerning the behavior of P(x) at the other stable flxed point, X4, it is either divergent 
(but integrable) for A;_/|/i(x4)| ^ 1, or it is continuous, 

(x) (fc+ + fc„)/i/:(x4)i-i 



lim P(x) 

a;\a;4 



lim P(x) 



L u{x,)[k./\fux,)\-iy 



(16) 



when k^/\f'_{xi)\ > 1. 

Imposing the normahzation condition for P{x), one obtains the corresponding asymptotic 
drift velocity: 



X4 



dx 



Xl 



X 



U{x)-f.{x) 



Uix)f.{x) 
k+ + k- 



dz sgn 



^3 



U{z)f-{z) 



X 



+ 



X exp ( - / dw{k+/Uiw) + k^/f4w)) ] + 



+ 



xi+L 



dx 



X4 



Uix)-f.ix) 



X 



k+ + 



(iz sgn 



X2 



h{z)-f-{z) 



+ 



f+{z) + f-{z) 



U{z)-f-{z) \2[U{z)-f4z)] 
X exp ( - [ dw{k+/f+{w) + k_/f_{w)) 



(17) 

\ J z / J 

which completes the discussion. 

All these results are suitable for a rapid generalization to cases when both dynamics have 
several pairs of fixed points. In particular, note that when there are several pairs of adjacent 
stable fixed points, the system is no longer ergodic. Depending on the initial conditions, the 
asymptotic motion of the particles is limited to one or another basin of attraction. These 
attractors are represented, on the real axis, by the intervals between such pairs of adjacent 
stable fixed points. We shall not go into further details here. 

IV. ROCKING RATCHET 



The general formulas that were derived above still involve triple integrals. One interesting 
case for which explicit results can be obtained is that of piece- wise constant functions /+ and 
To illustrate the results for this case, we will focus on one of the paradigms for Brownian 
motors, namely the rocking ratchet [2^. In this problem an overdamped particle is gliding 
in a deterministic sawtooth-like potential U (x) subject to an additional dichotomous forcing 
that switches between the values and — A_ . The equation of motion is of the form: 

A+ - A_ A+ + A^ 



-m 



(18) 
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FIG. 1: The potential U(x) used for numerical calculations on the rocking ratchet (/i=1.0, /2=0.5, 
L=1.0 and Li=0.25). 

with /(x) = —U'{x), and thus f±{x) = f{x)± A±. Without loss of generality, we assume 
that the mean value of the forcing is zero, i.e., A^/k^ = A^/k^, and, in order to fix ideas, 
we suppose that A^ > A^ > (the equality corresponds to a symmetric dichotomous noise). 
f{x) has a periodic "block- wave" profile: 



with f{x + L) = f{x), where L = Li + L2] also, /i and /2 are supposed to be two positive 
constants, with, e.g., /2 < fi- Although analytical calculation are done in full generality, 
all our numerical calculations were made for the potential represented in Fig. ^ i.e., we 
considered the case of a negative bias. In this case, the mean asymptotic velocity is deter- 
mined both by a ratchet-like effect (if present) and by the bias of the potential. Of course, 
one can also consider (as usually done when studying ratchet effect) the case of an untilted 
potential, that corresponds to the condition fiLi = /2i^2- 

We now proceed to a presentation of the results as they follow from the more general 
discussion presented earlier. Note that in the latter discussion f±{x) were assumed to be 
continuous. Therefore, in order to be able to use the results obtained previously, we shall 
first consider, for computational purposes, the following continuous, piece-wise linear force 




/i for X e [0, Li) 

/2 for X e [Li,Li + L2), 



(19) 
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profile: 



fix) 



fi for X G [0, Li - I) 

/i-(/i + /2)(x-Li + 0// for xe[L,~l,Li) 

-/2 for X G [Li, Li + L2 - /) 

-/2 + (/i + /2)(x-Li-L2 + 0A for xG [Li + L2-/,Li + L2). 



(20) 



Afterwards we sliall take the limit / [corresponding to the "block-wave" profile (jl9|) ] in 
the final expressions for the asymptotic probability density and the mean velocity. 



1. Strong forcing: no fixed points 



When A_ > /i and > /2, running solutions appear for both tilts and —A^, and 
there is no fixed point in any of the separate dynamics "+" and Applying the results 
obtained in Section IIIl A[ one finds the following expression for the probability density in 
the two subintervals of one period L: 



P{x) = < 



for X G [0,Li) 



jx) \ (/i + /2)A+A_(e^^-l) 1 



(21) 



for X G [Li, V) 



We used here the dimensionless quantities 

(A;+ + fc_)/iLi 



and 



(A;+ + fc_)/2L2 



(22) 



(A_-/i)(A+ + /i) ""^ (A_ + /2)(A+-/2)- 

Because of the first-order discontinuity of j{x) at x = Li, L, P{x) is also discontinuous at 
these points. 

The mean asymptotic velocity reads: 

-1 



{x) = L 



Li L2\ ^ A+A_ [ h + [l-e-^^ -e^^ + e^^- 



fi f2j k+ + k_ V fif 



(23) 



and it is determined by the interplay between the characteristics of the noise and those of 
the potential (in particular, its bias, if any). Some limit cases of interest include: 
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(i) When k± oo (but A± finite) one finds, of course, (x) = - the noiseless (determin- 
istic) result. 

(ii) The quenched noise-limit kj- ^ [i.e., a fraction A_/{A^ + A_) of the particles, 
chosen at random, are subjected to a constant external forcing A^, while the remaining 
ones are subjected to an external forcing —A^] results in the following mean velocity: 

^ fif2{Lj2 - L^h) - A+A_{LJ, - L2/2) - LfMA+ - A_) 



(x) 



(24) 



(L1/2 - L2/i)2 - L{A+ - A_)(Li/2 - L2/1) - • 

(iii) One obtains white shot noise ^ by taking the limit A_ ^ 00, A;_ — 00, such that 
A-/k- = A+/kj^ = A 7^ (y4+ and k+ being finite). The mean velocity is then 



{x)=L 



1 - e-*i - + e'f 



-\ -1 



(25) 



where 0i = fiLi/[X{A^ + fi)] and 02 = /2-^2/[A(A+ — / ?)]. In the case of a symmetric 
potential (/i = /2 and Li = L2), one recovers the result in |21|. 

Moreover, for very large values of A^ (and finite A) one recovers the noiseless limit, 
namely, (x) ~ (/i Li — /2 L2)/ L (that is strictly due to the bias of the potential). 

(iv) Finally, the white noise limit: A^ = A^ = A ^ 00 and k^ = k_ = k 00, with 
A'^/2k = D finite. Then 

— and 09 



JJ ^ D ' 

and the mean velocity ()23p takes the well-known form (see, e.g. |22|): 



(26) 



{x) = L 



fl f2j V /1/2 



g</'2-</'l — 1 



1 -1 



(27) 



In the case of an unbiased potential one has now 0i = 02 = 0, implying a zero mean velocity 
(as required by the second law of thermodynamics) and a probability density profile that 
assumes the Boltzmann form: 



P(x) 



/1/2 



(A + /2)D(1 - e- 



X < 



exp(— 0|a; — Lil/Li), for xe[0,Li] 



exp(— 0(x — Li)/L2), for xG[Li,L). 



(28) 



2. Intermediate forcing: two fixed points 



For intermediate forcing, there are two situations that might occur. 
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(1). When y4_ > /i but < /2, the "+" dynamics has two fixed points, namely, x = Li 
(stable, which corresponds to an asymmetric 6 peak in the probability density), and x = L 
(unstable, corresponds to a finite discontinuity in the probability density). One obtains from 
the results in Section ITllBI 

{x) 2A_(A+ + /i)(e<^i-l) 



P{x) 



L fc+/i(A+ + A. 
(x) 2A_(A+-/2)(e'^^ 
L k+f2{A+ + A. 



-5_(Li — x) 



-5+(x - Li] 



{xl 
Lfi 



{X 



A^-h 
-A_e'f"' 



exp[— — Lil/Li] — 1 



exp[-02(a; - Li)/L2] 



for X G [0, Li] 



for X G [Li, L). 



(29) 



A. + h 

Here 6±{x) are the half Dirac-delta functions [28]. A simple intuitive explanation for their 
appearance is related to the fact that in the limit of the "block- wave" force /(x), Eq. (fT^ . 
a particle can reach the stable fixed point x = Li (coming either from its left or from its 
right) in a finite time, and, once there, it will stay at this position for the reminder of the 
"+" dynamics. Thus, in time average, this will give rise to a finite weight at this precise 
point. The corresponding mean velocity is then: 

A_(A+ + /i)(e^^-l) 
/i 

A„(A+-/2)(e^^-l) 

/2 



+ 



+ 



1 



+ 



1 



K{A^ + A.) {K + k.)h 
1 1 



-1 



(30) 



Lfc+(A+ + A_) (A;+ + fc_)/2. 

The direction of the mean velocity is, as expected, the same as that of the "-" dynamics, 
i.e., (x) < in this case. 

One interesting limit in this case is that of white shot noise: A_ — * oo, A;_ ^ oo so 
that A_/k_ = A^/k^ = X finite. Our results are equivalent to those presented in il8| for 
the probability density (although written under a closed, compact form). One obtains the 
following expressions for the asymptotic velocity : 

-1 



-L 



'(L2 




A/2 


fj 



A(A+ + /i)2 {ef^ - 1) A(A+ - /2)2 {ef- - 1) 



A^fi 



(31) 
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0.5 1 1.5 

A 

+ 



FIG. 2: For the white shot noise Umit [eqs. (|25() . respectively H31() ]: The mean asymptotic velocity 
as a function of A-^- for various values of the transition rate A;+. The sawtooth potential is the one 
represented in Fig. ^ 

Figure 121 shows the variation of the average asymptotic velocity as a function of [according 
to Eq. ()25|1 . respectively (jS^)] various fixed values of the transition rate k+. Note the 
peak and the discontinuity in the slope of the velocity at = /2, connected with the 
appearance of the fixed points when decreaseas below /2. 

(2). When > f2 but A_ < fi, there are no fixed points in the "+" dynamics, 
but two fixed points in the "-" dynamics, namely x = Li (stable) and x = L (unstable); 
together with the discontinuous character of f{x) at these points, they lead, respectively, to 
an asymmetric (5-peak of P{x) at x = Li and to a first-order discontinuity of P{x) at x = L: 




+ 



{x) 2A+{A^~f\){e'^^-l) 
L k_MA+ + A_) 
(x) 2A+{A. + f,){e^--l) 
L k.UA+ + A_) 



5_(Li — x) 



6+{x - Li) 



exp[— — + 1 , for xG[0,Li) 



+ < 



(32) 



exp[-</)2(a; - Lx)lL'^ - 1 , for x G [Li, V). 



y Lf2 [A+ - /2 
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This leads to a mean velocity that reads: 



/i 

A+(A. + /2)(e^^-l) 

/2 



+ 



1 



+ 



1 1 

+ 



-1 



+ ' (fc+ + A;_)/2jj • ^^^^ 
Note that in this case (x) > 0, as dictated by the sign of the "+" dynamics. One can thus 
clearly realize the role of the ratchet effect if, e.g., one considers a potential with negative bias 
(like the one in Fig. Q): indeed, in this case the mean velocity is directed against the bias. 
The above expression simplifies for the case of a symmetric dichotomic noise Aj^ = A_ = A 
and = k_ = k, and reads: 



fi 



/2 



{A-hf{et>^-V^^{A + f, 



-1) 



2kfi 



2kfl 



(34) 



3. Weak forcing 



Finally, when A„ < /i and A^ < /2, there are no freely running solution in either of 
the dynamics. The points x = Li,L are fixed points for both "+" and "-" djTiamics. One 
concludes that the mean asymptotic velocity is zero, and all the particles are concentrated 
at the stable fixed point x = Li, i.e., P{x) = 6{x — Li). In this case additional thermal 
noise is needed to generate rectified motion. This problem has been solved for adiabatically 



slow forcing 10|, |ll | . 



^. Comparison with numerical simulations 

The Langevin equation (jlSp is simulated numerically using 100,000 particles and the 
probability distribution is averaged over 10 snapshots. The biased potential shown in Fig. 
H is used for all numerical calculations. Figures 01 and |3] show, respectively, the probability 
distribution P{x) and the average velocity (x) for the strong forcing case, and Figs. ElandlHl 
for the intermediate forcing case. Agreement between theory and simulations is very good. 
Discontinuities in P{x) aX x = Li and x = L are clearly seen in both cases. Figure El 
also illustrates the existence of a delta peak at x = Li. Note the phenomenon of current 
reversal when one decreases the amplitude of the perturbation (for a fixed k), i.e., when one 
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FIG. 3: Probability distribution P{x) for the strong forcing case {A^=A^=2.0, and k^=k_=5.0; 
see Fig. ^for the parameters of U{x).). The solid line shows the theoretical result 1)21(1 and the 
shaded histogram presents the results of numerical simulations. 




FIG. 4: Mean asymptotic velocity as a function of the transition rate k for the strong forcing case 
(see Figs. ^andlSlfor parameter values). The solid line shows the theoretical result (|23() and the 
open circles the results of numerical simulations. 
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X 



FIG. 5: Probability distribution P{x) for the intermediate forcing case {A^=A-=0.75, k-^-=k^=0.5; 
see Fig. ^for the parameters of U{x)). The sohd Une shows the theoretical result H29|) and the 
shaded histogram presents the results of numerical simulations. 

passes from strong forcing (no fixed points, cf. Fig. to intermediate forcing (fixed points 
appearing, cf. Fig. 



V. PERSPECTIVES 



The results obtained above reinforce the impression that dichotomous noise can be put 
on a par with Gaussian white noise as far as obtaining analytic results is concerned. Along 
this line of thought, we expect that one can obtain exact results for the first passage time 
moments when unstable fixed points are crossed. Also one can extend the calculations for 
the asymptotic diffusion coefficient in periodic system subject to additive Gaussian white 



noise 



23l 



2J| to the case of dichotomous forcing via its relationship to the first two moments 
of the first passage time Furthermore, the results presented here can be directly at 



plied to various other problems, including Stokes' drift and hypersensitive transport . 



suits pr 
ift y i 
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k 

FIG. 6: Particle velocity as a function of the transition rate k for the intermediate forcing case 
(see. Figs. nand|31for parameter values). The solid line is the theoretical result (|nn|) and the open 
circles result from numerical simulations. 

VI. APPENDIX 

From the point of view of the mathematical steps involved in solving Eqs. and ®, 
there are several cases that have to be considered separately: 

Case I, when the coefficient of p{x) in Eq. namely [f+{x) — f^{x)], has no zeroes in 
[0, L), one can solve this equation for p{x); 

Case II, when the coefficient of P{x) in i.e. [k^f^{x) + k^fj^{x)], has no zeroes in 
[0,L), one can solve ^ for P[x). Note that Cases I and II do not necessarily exclude one 
another. 

Case III, when both the coefficients of p{x) and of P{x) in Eq. Q have zeroes, hut not the 
same zeroes^ in [0,L), and finally 

Case IV, when f+{x) and f-{x) have (at least) a common zero in [0,L). In this case one 
can directly notice that there is no net flow of particles, i.e., J = 0, given the existence of 
the common flxed point that cannot be crossed in any of the two dynamics. This case will 
hence not be studied any further. 

A detailed analysis reveals that the results obtained in Cases I and II are completely 
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equivalent and apply also when [f+{x) — f-{x)] and/or [A;+/_(a;) + have different 

zeroes in [0, L], i.e., in Case III. We therefore focus exclusively on the main results pertaining 
to Case 1. 

In Case I one can solve Eq. (@)) with respect to p{x), i.e., 

J{k+ + k_) - [k+f_{x) + k.U{x)]P{x) 



p{x) 



(35) 



and by substitution in Eq. ^ one can obtain a closed first-order differential equation for 
P(x): 

[U{x)Ux)] P'{x) + - [f+{x)f^{x)] [In |/+(x) - f.{x)y; + 



+ \k^f-{x) + k_U{x)]} P{x) = J{{k+ + k-) + [f+{x) - /4a;)] 



[2[U{x)-f^{x)] 



(36) 



Recall that (...)' denotes derivation with respect to x. 

The crux of the problem resides in finding the correct solution to ()36p |26l |. A blind 
application of the standard method of variation of parameters leads to the familiar solution 

U(x)-f.(x) 



Pi.) = f 



[CG{x, xq) + K{x, xq] x)] 



(37) 



Uix)f.ix) 

where C is a constant of integration that arises from the general solution to the homogeneous 
part of Eq. (fHHjl . the second contribution is the particular solution of the full inhomogeneous 
equation , xq is an arbitrary point in [0, L), and we have defined the functions: 



G{u, v) 
K{u, v; w) 



exp 



dz 



+ 



k_ 



dz sgn 



Mz)-f_{z) 



k: + k. 



+ 



U{z) + f-{z) 



G{w,z) 



Mz)-Uz) \2[U{z)-f4z)], 

The point is now to get the correct integration constant C, or rather constants, since one 
can not apply the above solution at the points where the differential equation for P{x) is 
singular, which is precisely at the location of fixed points, i.e., at the zeroes of f±. 



A. No fixed points: /+(x)/_(x) 7^ in [0,L) 



In this case there is no fixed point in any of the alternating dynamics "+" and "-" . The 
usual procedure to determine C in Eq. can now be followed, i.e., we require periodicity 
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of P{x), i.e., P{x) = P{x + L), recalling that both f+{x) and /-(x) are periodic. This results 
in the expression Eq. ((Tj) of the main text. The mean velocity at the steady state, cf. Eq. 
(jSJ on the main text, follows from the normalization of P{x). 



B. One of the alternating dynamics has two fixed points in [0, L) 



The situation is entirely different, both physically and mathematically, when the system 
can cross unstable fixed points within the interval [0, L) in the long time limit. Consider 
the simple case described in Section Fill Bl in the main text. Clearly, the steady state results 
leading to Eqs. (jH) and (jHI) still apply, but the solution of Eq. © is more delicate than the 
"blind" integration that led to expression (jHTjl for P{x). Indeed, in the vicinity of the fixed 
points xi and X2, if we simply try to apply formula (jSZ|), we are led to the dependence 

P{x) ~ |x - xi,2r'"''+/t^+^"^''^l (38) 

For the case of the stable fixed point, P{x) is therefore continuous when > 
1, and divergent but integrable for ^ 1, a result that causes no conceptual 

difficulties. However, at the unstable fixed point X2, this procedure leads to an apparent 
non-integrable divergence, which is clearly unphysical, and mathematically improper in view 
of the requirement of normalization of P{x). 

As explained in detail in 19|, the fallacy lies in the assumption that a single integration 
constant C, see Eq. ^7\ . is valid throughout the region [0,L]. One solves the problem by 
choosing different integration constants in each of the separate intervals [0,xi), {xi,X2), and 
{x2,L] between the fixed points. There is exactly one choice of this constant valid for both 
{xi,X2) and {x2,L) such that the divergence at X2 is removed, namely C = —K{x2,xo;xo)] 
and another choice valid in the interval [0,xi) that ensures the required continuity and 
periodicity of P{x). The acceptable expression for the probability density is therefore found 
to be: 



P(x) 



{ M 

L 
L 



U{x)f.{x) 
f+{x) - f-{x) 



[K{L,X2;L)G{x,0) + K{x,0;xo)] for xe[0,xi) 



(39) 



K{X,X2]X) 



for X G (xi, L). 



/+(x)/_(x) 

These expressions can be further simplified if one takes as the basic period not 
[0,L], but + L]. Then the simple- looking, "compact" expression P{x) = 
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L 



U{x)-f4x) 



) [Eq. Q in the main text] holds throughout this new basic 

xi+L 



f+{x)f.{x 

period. With this choice of the basic interval, the normalization condition / P{x)dx = 1 
leads to the expression (jl2p for the mean velocity. 

P{x) as given above [Eq. ()39|) or Eq. Q] meets all the requirement enumerated in Sec. 
\n\ In particular, let us check its behavior at the fixed points xi and X2- In order to do this, 
it is useful to write Eq. (jH)) as: 

f+{x) - f-{x) 



Pi.) = M 



K{x,X2;0) , . 

|/+(x)|G(0,x)- ^ ^ 



For X = X2 (the unstable fixed point), K{x,X2; 0)/ G{0,x)} presents an indetermi- 

nacy of the type "0/0"; applying Hospital's rule one simply finds that P{x) is continuous in 
X = X2, and its value is given by Eq. (fT(H) . 

For X = xi (the stable fixed point), there are three situations that might appear, depend- 
ing on the value of k^/ |/:^(xi)|: 

(a) For /c+Z |/_^(xi) I > 1, both K{x,X2;0) and |/+(a;)| G{0,x) present a divergence for 
X \ xi, X y xi + L; therefore, P{x) presents an indeterminacy of the type "oo/oo". Ap- 
plying Hospital's rule one finds that P{x) has the same finite limit as x \ xi, x xi + L, 
as indicated in Eq. (fTTjl . 

(b) For /c+Z [/^(xi) I < 1, K[x,X2\^) is finite and nonzero at a; = X\,X\ + L, while 
|/+(x)| (7(0, x) as \x — X\\^^^ (respectively \x ~ X\ — V\~^^^\^^'^^^'^\\ There- 
fore, P{x) presents a power-law integrable divergence, P{x) ~ \x — x^^^^'^^^^^'^^^^'^ in the 
neighborhood of xi, respectively P(x) ~ \x — x\ — L\ in the neighborhood of 
xi + L. 

(c) Finally, for A;+/|/^(xi)| = 1, fC(x,X2;0) behaves like ln|x — xi| in the vicinity of x\ 
(respectively, like In |x — Xi — L\ in the vicinity of Xi + L), while |/+(x)| (^(0, x) has a finite 
limit in these points. Thus, P(x) has a "marginal", logarithmic-like divergence in these 
points. 

C. Each of the two alternating dynamics has two fixed points in [0, L) 

Suppose that each of the two dynamics, "+" and "-", has two fixed points (one stable 
and one unstable) in the interval [0,L), i.e., the situation described in Section Fill Bl Again, 
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if one blindly applies the result ()37|) for P{x), then in the vicinity of the fixed points 



P(x - xio) ~ la; - xio 



P{X - Xsa) 



-l-fc+/[/i(xi,2)] 



-l-k-/[f'_ix3,4)] 



(41) 



One again encounters the nonphysical non-integrable divergence at the unstable fixed points. 
The correct procedure is again to use the solution (|37|). but with different integration con- 
stants C in each of the open intervals between the fixed points. The latter constants are 
determined by the requirements (i) to remove the strong divergences at the unstable fixed 
points, by imposing the condition that the coefficient of the divergent term becomes zero 
at these points; (ii) to ensure continuity and periodicity of P{x); and finally, through the 
normalization condition, (iii) to determine the flow J and the mean asymptotic velocity (x) . 

Case a: ^ Xi (s) < X3 (u) < X4 (s) < X2 (u) < L 
The above "program" leads to the expression (fT^ in the main text for the probability 
density (in the appropriately chosen basic period + L]), and to the corresponding 

mean asymptotic velocity, Eq. (fTTj) . 

Case b: ^ Xi (s) < X2 (u) < X3 (s) < X4 (u) < L 
One gets the following expression for the probability density, again in the basic interval 
[xi,xi + L]: 

U{x)-f.{x) 



P(x) 



{ M 

L 

(i) 
L 



U{x)f.{x) 
U{x)-f_{x) 



K{x,X2,x),ioT X G {xi,Xs) 



K{x,X4,x),foT XE{x3,Xi + L) 



(42) 



/+(^)/-(^) 

with the good behaviors (continuity) at the unstable fixed points X2 and X4 and either 
continuity, or (integrable) divergences at the stable fixed points xi and X3; as well as the 
required periodicity for P{x). The mean asymptotic velocity is therefore given by: 



(x) 



dx 



U{x)-f.{x) 



/+(x)/-(^) 



K{x,X2,x) + 



x\+L 



dx 



^3 



/+(x)/_(s) 



K{ 



(43) 
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